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ABSTRACT 



Since its launch in 2004, the Swift satellite has monitored the X-ray afterglows of several hundred 
Gamma- Ray Bursts, and revealed that their X-ray light-curves are more complex than previously 
thought, exhibiting up to three power-law segments. Energy injection into the relativistic blast-wave 
energizing the burst ambient medium has been proposed most often to be the reason for the X-ray 
afterglow complexity. We examine 117 light-curve breaks of 98 Swift X-ray afterglows, selected for 
their high-quality monitoring and well-constrained flux decay rates. Thirty percent of afterglows have 
a break that can be an adiabatic jet-break, in the sense that there is one variant of the forward-shock 
emission from a coUimated outflow model that can account for both the pre- and post-break flux 
power-law decay indices, given the measured X-ray spectral slope. If allowance is made for a steady 
energy injection into the forward-shock, then another 56 percent of X-ray afterglows have a light- 
curve break that can be explained with a jet-break. The remaining 12 percent that are not jet-breaks, 
as well as the existence of two breaks in 19 afterglows (out of which only one can be a jet-break), 
suggest that some X-ray breaks arise from a sudden change in the rate at which energy is added to the 
blast-wave, and it may well be that a larger fraction of X-ray light-curve breaks are generated by that 
mechanism. The fractional increase in the shock energy, inferred from the energy injection required 
to account for the observed X-ray flux decays, may be anticorrelated with the GRB prompt output, 
whether the X-ray break is a jet-break or an energy-injection break. That anticorrelation can also 
be seen as bursts with a higher energy output being followed by faster fading X-ray afterglows. To 
test the above two mechanisms for afterglow light-curve breaks, we derive comprehensive analytical 
results for the dynamics of outflows undergoing energy injection and for their light-curves, including 
closure relations for inverse-Compton afterglows and for the emission from spreading jets interacting 
with an wind-like ambient medium. 
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1 INTRODUCTION 



It was also realized early on that there should be a sub- 
stantial diversity in the steepness of the afterglow flux decay, 
even when the spectral slope is constant, if the fundamen- 
tal forward-shock parameters vary as the shock traverses 
the ambient medium. Rees & Meszaros (1998) derived the 
expected power-law flux decay for the forward-shock syn- 
chrotron flux in the case of an outflow whose kinetic energy- 
per-solid angle dE/d£l is not constant in time, either because 
the shock's energy increases due to fresh ejecta arriving at 
the blast- wave or because dE/dO, is anisotropic. Two impor- 
tant microphysical parameters (ec and ee), which quantify 
the fraction of the post-shock energy acquired by electrons 
and magnetic fleld, also determine the intensity the after- 
glow synchrotron flux, and their variability could also yield 
a variety of decay indices a. 



The power-law temporal decays displayed by the radio, op- 
tical, and X-ray fluxes of GRB afterglows [F^ oc t~°') have 
provided strong support for the prediction (Meszaros & Rees 
1997) that afterglows arise from the relativistic blast-wave 
generated as the GRB outflow ejecta interacts with the am- 
bient medium. If the GRB outflow is a coUimated jet then 
a break in the smooth power-law decay of the afterglow is 
expected for both either a conical jet (Panaitescu, Meszaros 
& Rees 1998) or a sideways-spreading jet (by Rhoads 1999) . 
The breaks found in the flux decay light curves of most of 
the well-sampled afterglows that were localized by the Bep- 
poSAX satellite seemed to provide solid observational sup- 




et al 2001). 
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Swift measurements of GRB X-ray afterglows pro- 
vide support for that shock parameters vary. The standard 
forward-shock model with constant energy and shock micro- 
parameters predicts a closure relation between the tempo- 
ral decay index a and the slope jS of the X-ray continuum 
{F„ oc f of the form a = 1.5/3 -|- const. The early X-ray 
flux of some Swift afterglows decays slower than expected 
for the measured spectral slope fi, which was attributed to 
a variable kinetic energy of the forward shock (Nousek et al 
2006, Panaitescu ct al 2006, Zhang ct al 2006). Moreover, 
the flux decay indices a displayed by Swift X-ray afterglows 
after the slow decay phase are not correlated with the X-ray 
spectral slope /3 (Panaitescu 2007), as would be expected 
from the a — P closure relation for the forward-shock emis- 
sion. 

When optical light-curves are taken into account, the 
standard jet model fares even worse, being inconsistent with 
practically all Swift afterglows (Liang et al 2008). While 
more than half of the well-sampled optical afterglows dis- 
play coupled optical and X-ray light-curves (Panaitescu & 
Vestrand 2011), exhibiting similar flux decay indices and/or 
simultaneous light-curve breaks, many other afterglows show 
a decoupling of their optical and X-ray light-curves, with a 
light-curve break occurring only in the X-ray. At least four 
possibilities for this light-curve decoupling, which is not pre- 
dicted by the standard jet model, have been proposed: i) the 
X-ray afterglow emission arises occasionally from the same 
not-yet-understood mechanism which produces the prompt 
burst emission (Ghisellini et al 2007), ii) there is a dominant 
contribution of the reverse-shock to the afterglow emission 
(Uhm & Beloborodov 2007) , in) the forward-shock emission 
is reprocessed by bulk-scattering and/or inverse-Compton 
scattering outside the forward-shock (Panaitescu 2008), or 
iv) the optical and X-ray afterglow flux arise from different 
parts of the forward-shock, from different jets (e.g. Racusin 
et al 2008). 

In this work, we assume that Swift X-ray afterglows 
arise from the forward-shock, and attribute the failure of the 

standard jet model in accounting for the X-ray flux decay 
indices to the Eissumption that the forward shock energy is 
constant. 

Evolving micro-physical parameters (e.g. loka et al 
2006) could, in principle, explain the measured X-ray flux 
decay indices, but the existence of light-curve breaks re- 
quires an unnatural, sudden change in their evolution. For 
that reason, we maintain in this work the basic assumption 
of the standard jet model that micro-physical parameters do 
not evolve. 

We also assume that the ambient medium with which 
the post-GRB outflow interacts is uniform, without discon- 
tinuities that could yield breaks in the afterglow light-curve 
at a frequency. A notable model thus excluded is that where 
the afterglow-producing outflow interacts with a prior out- 
flow (loka et al 2006) whose density increases radially out- 
ward up to its termination shock, yielding a slower decay of 
the afterglow flux than expected for a homogeneous ambient 
medium, followed by a break to a steeper decay. 

One reason for relaxing the constant shock energy con- 
straint is the observational evidence for X-ray flares in Swift 
afterglows (e.g. Chincarini et al 2007). The short timescale 



of those flares indicates that the GRB progenitor produces 
relativistic ejecta after the burst phase, whose arrival at the 
decelerating forward-shock could increase substantially the 
blast-wave's energy. Our hypothesis is that such an injec- 
tion of energy in the forward-shock is, in general, a more 
persistent and continuous process. However, we note that 
the fluence of X-ray flares is about 10 percent that of the 
burst (Falcone et al 2007), hence the late outflow powering 
the flares could carry a dynamically-important energy only 
if that outflow is much less efficient than the GRB outflow 
in producing X-ray emission. 

We derive in §2 the closure relations a — /3 for a con- 
tinuous energy injection in the forward-shock. For the ease 
of extracting from observations the required injected power 
history, we assume that the forward-shock energy increases 
as a power with observer-time, E oc t'', which is motivated 
by that the resulting afterglow flux should be power-law in 
observer time. This assumption circumvents a calculation of 
the kinematics for the ejecta - forward-shock interaction. 

In §3, the closure relations for the forward-shock emis- 
sion with energy injection are applied to a set of Swift af- 
terglows with well sampled light-curves. There, we exam- 
ine two plausible scenarios for the origin of afterglow X- 

ray light-curve breaks (break in energy injection, jet- break), 
three dynamical regimes of the forward-shock (spherical ex- 
pansion, conical jet, spreading jet), and two radiation pro- 
cesses that may yield the afterglow X-ray flux (synchrotron, 
inverse-Compton) . 



2 ANALYTICAL FLUX DECAYS FOR 
FORWARD-SHOCK MODELS WITH 
ENERGY INJECTION 

There are three mechanisms that could produce a change 
in the flux power-law decay index a, i.e. an afterglow light- 
curve break. 

(1) Passage of a spectral break. X-ray light-curve breaks 

are, generally, not accompanied by a spectral evolution 
(Nousek et al 2006; Liang, Zhang & Zhang 2007; WiUingale 
et al 2007), hence their origin in the passage of the spectral 
feature is largely ruled out, but there axe other reasons why 
such an X-ray light-curve break origin is unlikely. The slower 
evolution of the cooling frequency Vc (the synchrotron char- 
acteristic frequency at which radiate the electrons that cool 
radiatively on the hydrodynamical timescale) leads to ex- 
tended X-ray light-curve breaks when Vc crosses the 0.3-10 
keV band. The peak frequency z/i of the vFu synchrotron 
spectrum decreases faster, but it is likely to be lower than 
the X-rays even at early times, hence it cannot cross the 
X-ray band during afterglow observations. 

(2) Jet-break. A light-curve break occurs at the 'jet- 
break time', when the outflow Lorentz factor F drops below 
the inverse of the jet half-angle dj. For a 'conical jet' that 
does not spread laterally (as in the case of an outflow core 
embedded in an envelope that prevents the spreading of the 
core but yields a negligible emission) , the steepening of the 
X-ray flux decay is solely due to that, after the jet-break 
time, the average kinetic energy per solid angle dE/dil over 
the visible cone of angle 6 = T~^ (where relativistic Doppler 
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boost is effective) is decreasing with observer-frame time. 
For a 'spreading jet' (as expected for a top-hat distribution 
of dE/dQ), in addition to the above-described geometrical- 
relativistic effect, there is a second contribution to the X-ray 
flux decay steepening from the enhanced jet deceleration 
produced by the jet lateral spreading. 

(3) Energy injection break. If energy injection occurs 
routinely in GRB afterglows, then a light-curve break could 
be produced by a change in the injected power. The dynam- 
ical behaviours of the forward-shock model to be consid- 
ered are: i) a 'spherical' outflow, i.e. a wide jet whose edge 
becomes visible to the observer after the last X-ray mea- 
surement (at tend), which occurs if jet half-angle satisfies 
0j > l/r(tend), and ii) a narrow jet whose edge is already 
visible at the first X-ray observation (at tstart), occurring 
when Oj < l/r(tstart). For the latter model, we will con- 
sider both a conical jet and a spreading jet. These dynam- 
ical models are applicable also to X-ray afterglows without 
light-curve breaks, in which case the energy injection law 
must be unevolving (a single index e). 

We will investigate the ability of the energy-injection 
break and jet-break models above to accommodate the ob- 
served X-ray flux decays before and after the break, with 
allowance for all three dynamical models (spherical out- 
flow, conical jet, spreading jet), for two types of ambient 
medium (homogeneous and wind- like), two radiation mech- 
anisms dominating the X-ray flux (synchrotron and inverse- 
Compton), and two possible locations of the X-ray relative 
to the cooling frequency i>c but, first, a word about their 
rationale. 

Wind-like media, having a n ~ 0.1 - l(r/10^®cm)-2 
particle density, are expected at 0.1-1 pc from the source 
of long GRBs, i.e. where the 100 s — 1 Ms afterglow emis- 
sion is produced, because the progenitors of such bursts are 
massive Wolf-Rayet stars that drive powerful winds. The 
termination-shock of such winds is of several pes (Garcia- 
Scgura ct al 1996), thus a homogeneous medium at the 
afterglow location seems unjustified. However, if the GRB 
progenitor resides in a highly-pressurized bubble, the wind 
termination-shock radius could be as low as 0.5 pc (Cheva- 
lier et al 2004), with a further reduction of that radius if 
the WR wind is weak, as expected along the rotation axis 
(Ramirez- Ruiz et al 2005), or if the WR star moves suffi- 
ciently fast relative to the ISM (van Marie et al 2006), thus 
the afterglow could occur in shocked wind, which is quasi- 
homogeneous. 

Synchrotron emission is generally believed to produce 
all the afterglow emission, but the likelihood of inverse- 
Compton emission being dominant is higher at higher pho- 
ton frequency, like X-rays. Admittedly, that requires a high 
environment density, perhaps ten times higher (Panaitescu 
2011) than implied by the mass- loss rates and terminal wind 
velocities measured for Galactic Wolf-Rayet winds (Nugis & 
Lamers 2000). 



i^Fv spectrum peak), and fc (cooling frequency), and from 
the shape of the spectrum at frequencies higher than its 
min(i/i, i^c) peak: F^, oc i/'^^^ for fc < < vi, F^ oc 
for Vi < V < and F,^ oc for raax.{vi, Vc) < v, where 

p is the power-law index of the electron distribution with 
energy produced by the shock (dN/de oc e^^). Therefore, 
the flux received at observing frequency v satisfies F„ cx 
F^vl"-^'"^ for Vi < u < V, and oc F^v^'^'^'^ vl'^ for 
max(fi,i^c) < V [v^ < V < Vi implies an X-ray spectrum 
with /3x = 1/2 that is harder than measured for nearly all 
Swift X-ray afterglows). 

For a spherical outflow (or a jet before the jet-break 
time tj) and for synchrotron emission, _Fp cx BMV, Vi oc 
7fBr, Vc oc ■y^BF. Here M is the mass of the ambient 
medium swept-up by the forward-shock (i.e. the number of 
radiating electrons); B oc (eBUj^^^T is the post-shock mag- 
netic field, whose energy density B^ /{8tt) is assumed to be a 
fraction en of the post-shock energy density u' = AF^nmpC^ 
(from the shock jump conditions); n is the external density 
at the shock's location; 71 oc eeF is the typical electron ran- 
dom Lorentz factor, parametrized as a fraction to of the pro- 
ton random Lorentz factor jp ~ F; and 7c oc [B'^t' {Y + 
is the random Lorentz factor of the electrons whose radia- 
tive cooling timescale equals the comoving-frame shock's age 
t' ~ f dr/{cT), r being the outflow radius and Y the Comp- 
ton parameter. We assume the more likely situations of ei- 
ther electrons cooling mostly through synchrotron emission 
(Y < 1) or, if inverse- Compton is dominant (Y > 1), that 
z/c < I'i, in which case Y ~ (ee/ee)'^'^^ is constant and does 
not affect the evolution of Vc- 

For a jet (conical or spreading) after the jet-break time 
and for synchrotron emission, we have Fp oc BMF^, with Vi 
and Vc as above. Thus 



u^'""'^ oc v/^F" 



{^MV^ it<tj) 
{t>tj) 



v'f'"'^ oc (n^'^V't 



3/2-p2 ,/2n -1 



(1) 

(2) 

(3) 



The characteristics of the inverse-Compton spectrum 
are immediately related with those of the synchrotron spec- 
trum: ~ TeiF^'"''"^ (where Tei oc M/r^ is the opti- 
cal thickness to electron scattering of the shocked ambient 
medium), i/'''^'' ~ 7?^'i'^''°', and Vc ' — thus 



F('«) oc 



K^'C' oc 0^r« 



2.2 Outflow dynamics 



(4) 

(5) 
(6) 



2.1 Spectral characteristics 

The a—/? closure relations expected for a given variant of the 
forward-shock model are derived from the evolution of the 
spectral characteristics Fp (peak flux), vi (frequency of the 



The above equations for the spectral characteristics show 
that their evolution and, thus, the resulting flux decay in- 
dex a, is set by the forward-shock dynamics: T(t), M{t), and 
r{t). The energy-per-particle behind that shock is e/n = F, 
hence the lab-frame energy of the forward-shock is F^Mc^ 
and its dynamics is given by F^M = E/c^ = (MoTo){t/toY , 
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where MoFq is the kinetic energy of the cjccta (all trans- 
ferred to the forward-shock) at the time to when the as- 
sumed power-law energy injection {E oc t'^) begins. The 
other equations for the forward-shock dynamics are dM = 
Ar^"{r9)^dr, with p — Ar^" the density of the ambient 
medium at radius r (s = for a homogeneous medium, 
s = 2 for a wind) and dO = Csdt' /r = dr/{\/irT) the in- 
crease of the jet opening at the comoving-frame sound speed 
(cs = c/a/3), and dt = dr/{cV'^) for the observer-frame ar- 
rival time t of photons emitted along the direction toward 
the observer. 

Defining / = M/Mq, x = r/r with f = 
(GMoFo/ttA)^'''''^^' (the radius where the lateral spreading 
almost doubles the initial opening of an adiabatic jet in- 
teracting with a homogeneous medium), and t = t/to, the 
above equations for the outflow dynamics become: 



To e 



da; 

de 

dx 

dr f 
dx ^ 2ctor2 



VSxT 



(7) 
(8) 
(9) 
(10) 



Equation (10) allows the conversion of the jet dynamics 
[r(r), M(r), d{r)] from lab frame to observer frame [T{t), 

M{t), eit)]. 

For a spherical outflow, a conical jet, or a spreading 
jet before the jet-break time, 9 is constant and the above 
equations can be solved easily, leading to 



r oc r 



e)/(8-2s) 



A^^^(3-.)(i+e)/(4-.)(^ nr^) 



r oc t 



(H-e)/(4-3) 



^'^^(5-.+e)/(8-2.)(^ r/r) 



(11) 
(12) 
(13) 
(14) 



which can be replaced in the equations for the synchrotron 
and inverse-Compton spectral characteristics (equations 1 3 
and 4-6) to obtain their evolutions and, then, the power-law 
index of the afterglow flux decays given in lines 1, 2, 4, and 
5 of Table 1 and Table 2, respectively. 



2.3 Spreading jets 

2.3.1 Adiabatic jet, homogeneous medium 

For a jet spreading laterally, the dynamics equations can 
be solved analytically in the case of no energy injection 
and for a homogeneous medium. That has been done by 
Rhoads (1999), who obtained an exponential jet deceleration 
for a lateral spreading given by 6* = c^t' /r: T oc exp(— fcr), 
6 oc cxp(fcr)/'r, / oc cxp(2fcr) at r > k^^ , where k — 1/f. Be- 
cause the differential lateral spreading is "on top" of a coni- 
cal jet and not "over" a cylindrical jet (as in Rhoads' formal- 
ism), the correct lateral expansion is dO = c^dt' jr. For e = 



and s = 0, equations (7)-(10) (which use the latter prescrip- 
tion for jet spreading) have the solution F oc cxp(— fcr)/r, 
Q oc exp(fcr), / oc exp(2fcr) for r > f, which can be 
shown to lead to the same flux power-law decay indices as for 
Rhoads' prescription because the increase of the jet radius 
is only logarithmic in observer time (note that, for a nearly 
constant jet radius r, dB = Cgdt' /r implies 6 = Cat' /r). 

The above exponential solutions for the dynamics of 
an adiabatic jet interacting with a homogeneous medium 
can be expressed in observer time by integrating t = t{f) + 
/J' dr/(2cF^), which leads to t ~ f/(4cF^) for r > f. Thus 



1/2 



from where it can be shown that 

M = ^ oc E^/'n'^h 
rc 



, 2 /I. \^/^ /E\ 



1/2 



/2 



(?) «itJ 

where E is the jet energy and f = (GE/nnmpC^y^^ . 



(15) 

(16) 
(17) 
(18) 



2.3.2 Energized jets 

In the more general case of energy injection and/or for a 
wind-like medium, the solution to the dynamics equations 
(7)-(10) cannot be found analytically. They can be cast as 
logarithmic derivatives 



din/ r 



2-s e+1 



dln6» r 

oc 



l-l-e/2 



dlnr 



d In X r*^ 

oc 



d In r 



fx 



(19) 
(20) 
(21) 



and integrated numerically, which showed that the left-hand 
sides are (asymptotically) constant, i.e. the dynamics [M{t), 

0(t), r(t)] of a spreading jet undergoing energy injection ac- 
cording to i5 oc f*^ is a power-law in observer time. That 
dynamics is between the extreme cases of an adiabatic jet 
interacting with a homogeneous medium (which yields the 
fastest deceleration) and a spherical outflow with energy in- 
jection (which yields the slowest deceleration). Given that, 
in both of those two extreme cases, the dynamical quanti- 
ties have a power-law evolution with time, it is not surprising 
that an intermediate case manifests that feature too. 

Replacing power-law solutions for f , 6, and x in equa- 
tions (19)-(21), allows the calculation of the temporal ex- 
ponents by requiring that the right-hand sides are also con- 
stants, i.e. without a r dependence. The result is 



Moct" 



■(oc nr^e^) 



oc t 



5(1-3^) 



(22) 
(23) 
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Table 1. Evolution of the spectral characteristics (Fp = spectral peals flux, f i = energy at which radiate the electrons with typical post- 
shock energy, Vc = cooling frequency) with observer time t for the (synchrotron and inverse-Compton spectra, for two types of ambient 
media (homogeneous and wind), for three models for the forward-shock dynamics ('spherical' = wide jet, with a jet-brealf after last 
observation, 'conical jet' = confined jet, observed after jet-break, 'spreading jet' = laterally spreading jet, observed after jet-break) , and 
a forward-shock energy increasing as E oct". 



Model Homogeneous medium Wind-like medium 
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syn 


- spherical 


e 


4(3- 


e) 


4(1 + ^) 


4(1- 


e) 


4(3 


-e) 


^(1 + e) 


syn 


- conical jet 


i(5e-3) 


4(3- 


e) 


-i(l + e) 


-(1- 


e) 


4(3 


-e) 


i(l + e) 


syn 


- spreading jet 


|e-l 


-1(3- 


e) 


-h 


-(1- 


e) 


-(2 


-e) 


e 


iC - 


spherical 


3(5e-H) 


-1(3- 


-e) 


4(5e + l) 


-1 




-(2. 


-e) 


2 + e 


iC - 


conical jet 


i(3e-l) 


-1(3- 


-e) 


4(5e + l) 


4(3- 




-(2 


-e) 


2 + e 


iC - 


spreading jet 


h 


-(3- 


e) 


-(fe-1) 


— e 




-(3- 


-2e) 


l + 2e 



Table 2. Afterglow flux power-law decay index a {Fi, oc t~°') for the models listed in Table 1. p is the power-law index of the electron 
distribution with energy behind the forward-shock; dN/dt oc e~^, determined from the spectral slope /S: p = 2f3 + 1 for vi < u < Vc and 
p = 2/3 for i/; , z/c < t". 



Model Homogeneous medium Wind-like medium 
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r oc (24) 

It is worth mentioning that the above dynamics of spread- 
ing jets undergoing energy injection can be recovered from 
the dynamics of a spreading adiabatic jet interacting with 
a homogeneous medium (equations 16, 17) if one replaces 
E cc t"^ (the prescription for energy injection) and r oc T^t. 
The latter relation is trivial for an outflow decelerating as a 
power-law with radius (as for a spherical outflow or a conical 
jet) and is approximately satisfied for an exponential decel- 
eration (as for an adiabatic spreading jet interacting with a 
homogeneous medium), which suggests that it should also 
stand in any other case of an intermediate-strength deceler- 
ation (that is verified with the numerical integration of the 
jet dynamics). Either from equations (22) and (24)or by re- 
placing E and n oc r~" in equations (15) and (18) with the 
above prescriptions, it can be shown that 

roc<"5(i-3^^> (25) 

octi''+3^'(ocr/r) (26) 
for a spreading jet and with energy injection. 



By substituting equations (22)-(26) in the relations for 
the synchrotron and inverse-Compton spectral characteris- 
tics (equations 1-3 and 4-6), one can derive the evolutions 
from lines 3 and 6 of Table 1 and the power-law indices of 
the afterglow fiux decay from lines 3 and 6 of Table 2. 

We acknowledge that the dynamics of spreading jets 
undergoing energy injection was derived using a simplistic 

toy model for the jet lateral spreading, where the jet was 
assumed to have a uniform kinetic energy per solid angle 
dE/dil,, although both lateral spreading and energy injec- 
tion could/should introduce a substantial angular gradient 
of dE/dO,, decreasing with angle measured from the jet's 
spine. 



3 ENERGY INJECTION IN SWIFT X-RAY 
AFTERGLOWS 

3.1 Sample selection 

Out of the 634 afterglow light-curves listed at the Swift- 
XRT light-curve repository (Evans et al 2007) posted until 
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Table 3. Light-curve breaJcs of well-sampled Swift X-ray afterglows. GRBs in bold-face have afterglow light-curve breaks that can be 
accounted for by an adiabatic jet-break (without energy injection). The GRBs listed in italics have a light-curve break that cannot be 
explained with a jet-break uiidorgoinji an oiiorgy iDjoctioij continuous across the liioak 
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Table 3. (continued) 
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(1) : X-ray spectral slope (90 percent CL) 

(2) : epoch of first Swift X-ray measurement 

(3) : index of X-ray flux power- law decay at tstart — tb (90 percent CL) 

(4) : epoch of X-ray hght-curve break 

(5) : index of X-ray flux power-law decay at tt ~ tend (90 percent CL) 
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the end of year 2011, we have retained only those for which i) 
a light-curve break is obvious or seems possible, ii) a power- 
law light-curve segment is monitored for longer than a factor 
5 in time before and after the break, and Hi) the spectral 
slope /3x listed at the Swift-XRT spectra repository (Evans 
et al 2007) has a 90 percent confidence level (CL) less than 
0.16. 

The hght-curvos have boon pioco-wiso fit with power- 
laws, using measurements sufficiently far from the break 
time, to capture the asymptotic power-law flux decay. Shal- 
low breaks for which the data could be fit satisfactorily (i.e. 
reduced — 1) with a single power-law were excluded. The 
90 percent confidence level {a(a^)) was determined from the 
variation of around the best power-law fit. 

As shown in Table 2, for an adiabatic outflow (e = 

0) , the power-law index is Qx = 3p/4 + const = 3/3x/2 + 
const, with the coefficient of /3x being as large £is 9/8 for 
other models. Thus, an equal importance of measurement 
uncertainties a(Qx) and <t(/3x) < 0.16 requires that we retain 
only afterglows whose decay indices can be determined with 
a 90 percent CL smaller than 1.5cr(/3x) — 0.24. 

The above two basic selection criteria led to the sample 
of 98 afterglows and 117 X-ray light-curve breaks listed in 
Table 3. Forty-four afterglows have temporal coverage start- 
ing from the burst end (for them, we determine the fractional 
increase of the outflow energy during the afterglow phase) 
and a measured redshift, which allows the calculation of the 
GRB output and of the afterglow jet energy. 

The light-curves of the afterglows listed in Table 3 are 

morphologically diverse. A minority display a continuous, 
slowly-decaying power-law flux from the burst phase to the 
afterglow, the afterglow beginning being chosen when the 
flux fluctuations disappear. The majority of GRBs display a 
sharp decay after the prompt phase (the GRB tail) followed 
by a much slower decay and spectral hardening; the after- 
glow beginning is chosen at the epoch when the slower flux 
decay starts. 

The closure relations a — /3 listed in Table 2 are used 
to calculate the required energy injection indices from the 
measured pre- and post-break X-ray flux decay indices a and 
the X-ray slope /3 (which gives the electron index p) . For a 
model to be validated, the resulting energy injection index e 

1) should be positive (i.e. e + (Te > 0, with the uncertainty as 
calculated from the propagation of cr(ax) and a(/3x)), and ii) 
for jet-break models, should lead to a decelerating outflow, 
so that the jet boundary can become visible to the observer 
at some time. Prom equation 11, this condition is e < 3 — s. 



3.2 Adiabatic jets 

For adiabatic jets (e = 0), the obscrvablcs a^i and ax2 yield 
two constraints. As there are no free model parameters, aside 
the various possibilities related to the jet dynamics, location 
of the cooling frequency, ambient medium type, and domi- 
nant X-ray emission process, the model is over-constrained. 
Thus, it should not be surprising that only a small fraction 
of Swift X-ray breaks satisfy |e| < ae both before and after 
the break, i.e. can be explained by any single variant of the 
adiabatic jet break model (upper part of Table 4). 



Table 4. Number of breaks that can be accommodated by an 
adiabatic jet-break, a jet-break with steady energy injection, or 
by an energy-injection break, out of the 117 breaks listed in Table 
3. Models are identified by the origin of the X-ray light-curve 

break ("J" = jet-break, "EI" = energy-injection break), radiation 
process dominant in the X-rays ("syii" = synchrotron, "iC" = 
inverse-Compton) and by the jet dynamics ( "cjct" = conical jet, 
"sjet" = spreading jet, "sph" = spherical - see Table 2). 



MODEL HOMOG MEDIUM WiND-LIKE MEDIUM 

U < Uc Uc < U U < Uc Uc < U 



adiabatic jet-break 



J-syn-cjet 


3 


11 





9 


J-syn-sjet 





7 





7 


J-iC-cjet 


2 


15 





9 


J-iC-sjet 


2 


7 





7 


jet-break with 


or without 


steady EI 




J-syn-cjet 


48 


47 


12 


29 


J-syn-sjet 


19 


39 


15 


42 


J-iC— cjet 


64 


52 





31 


J-iC— sjet 


67 


18 





16 


energy-injection break 


El-syn-cjet 


111 


98 


113 


77 


El-syn-sjet 


115 


106 


115 


105 


EI— syn— sph 


77 


31 


106 


31 


El-iC-cjet 


114 


4 


116 


80 


El-iC-sjet 


114 


47 


114 


54 


El-iC-sph 


105 


37 


115 


36 



Evidently, Swift X-ray afterglows may not be all of the 
same forward-shock model type. Even when we take into 
account all these possible variants, we find that adiabatic 
jet-break model can explain only 30 afterglows (indicated in 
boldface in Table 3), i.e. about 30 percent of the sample con- 
sidered here. Six of those 30 breaks occur in afterglows with 
two light-curve breaks, the break compatible with an adia^ 
batic jet-break being always the second. For all six, there is 
always a spherical energy-injection break model with e = 
after the break that accounts for the first break, and com- 
patible with the jet-break model, i.e. the two break models 
have the same location of cooling frequency relative to the 
X-rays, the same type of external medium, and the same 
X-ray emission process. 



3.3 Energized jets 

3.3.1 Jet-breaks 

For the jet-break model, we require that the X-ray light- 
curve break arises solely from the outflow coUimation, with- 
out any contribution from a break in the injected power, i.e. 
we require a constant energy injection index e across the X- 
ray light-curve break: (61 — 62)^ < a^^+a^^- Thus, the model 
is still over-constrained, with two observables (oxi and 0x2) 
that must be accommodated by a single model parameter 
(e). Considering that the jet-break model with steady energy 
injection injection is less over-constrained than the adiabatic 
jet model, it may not be surprising that the former model 
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logioE^(erg) 

Figure 1. Fractional increase of the jet energy from (source 
frame) 100 s to 100 ks (-Eend /^start ) versus GRB output at 10 
keV — 10 MeV (E^) and whose X-ray light-curve breaks can be 
accommodated by a conical jet-break with steady energy injection 
(other model details: X-ray higher than the synchrotron cooling 
frequency, external medium is homogeneous). Black points are for 
a model where synchrotron is the dominant X-ray radiative pro- 
cess, red points for an inverse-Compton model. These are the af- 
terglows listed in Table 3 for which a redshift is known and whose 
light-curve breaks (first break, if two exist) can be accommodated 
by the models described above. In natural scale, Eend/^'start and 
£'end/£'start are Uncertain by up to 50 percent, hence their loga- 
rithms are uncertain by up to 0.20 . 



can account for a larger number of Swift breaks (mid part 
of Table 4). 

With allowance for all jet-break model variants listed in 
Table 4 (i.e. for either location of the cooling frequency rela- 
tive to the X-ray, for either radiation process, for either type 
of medium, and either post jet-break dynamics, conical or 
spreading), the jet- break model with steady energy injection 
(including adiabatic jet-breaks, for which the energy injec- 
tion parameter e is consistent with zero) can account for 103 
of the 117 X-ray light-curve breaks considered here (the un- 
explained breaks are shown in italics in Tables 3). 37 of those 
potential jet-breaks occur in afterglows with two light-curve 
breaks; with the exception of the first break of afterglow 
081222, all other 36 jet-breaks can be paired consistently 
with an energy-injection break for the remaining break, i.e. 
the two break models have the same medium stratification, 
the same X-ray radiative process, the same location of X- 
ray relative to the cooling frequency, and compatible ejecta 
dynamics. 

Because only one break can arise from the outflow coUi- 
mation and 19 afterglows in our sample have two light-curve 
breaks, a more correct way of assessing the success of the 
jet-break model with steady energy injection is to say that 
it can explain at least one of the breaks observed for 85 af- 
terglows out of 98 or, conversely, that it cannot account for 
any of the X-ray light-curve breaks of 12 percent of after- 
glows. Excluding the 30 breaks (of 30 afterglows) that can 
be accounted for with an adiabatic jet-break, the jet-break 



model with non-zero energy injection can explain one X- 
ray light-curve break for 55 of the 98 afterglows considered 
here, which argues strongly that energy injection may be 
often present in Swift afterglows. 



3.3.2 Energy injection breaks 

Given that 12 percent of afterglows have an X-ray light- 
curve break that cannot be explained with a jet-break even 
when a steady energy injection is allowed, and that some 
afterglows display two breaks, out of which only one can 
be a jet-break, it is natural to speculate that some X-ray 
light-curve breaks must arise from a change in the power at 
which energy is added to the forward-shock. In the extreme, 
one could speculate that all X-ray breaks could be such en- 
ergy injection breaks, but the large fraction of breaks that 
can be explained as jet-breaks argues against that extreme 
speculation. 

The broken power-law energy injection model is not 
over-constrained, as there are two observables (ctxi and 0x2) 
for two model parameters (ei and 62), the only restriction 
being that e > 0. Thus, it should not be surprising that 
there are variants of the energy-injection break model that 
can account for most of the Swift X-ray breaks (lower part 
of Table 4). There is only one X-ray break that cannot be 
explained with an energy injection break: GRB afterglow 
060607. From Table 2, it can be seen that, for /3x < 3/2, the 
steepest possible decay for an adiabatic outflow is that for 
the inverse-Compton emission from a conical jet interacting 
with a wind, and for < fc: Ox = p + 1/2 = 2/3x + 3/2 
*. For GRB afterglow 060607, /3x = 0.61 ± 0.06, hence the 
steepest possible decay for the forward-shock emission has 
Qx = 2.72 ± 0.12, which is slower than the measured post- 
break 0x2 = 3.16 ± 0.07. 



4 IMPLICATIONS OF ENERGY INJECTION 

4.1 Afterglow decay rate — GRB output 
correlation 

After having identifled the power-law energy injection that 
accommodates each of the X-ray light-curve breaks, one can 
calculate the increase in the forward-shock energy from flrst 
to last X-ray measurement. The interesting result is that 
the fractional increase of the forward-shock energy over a 
source-frame time-interval common to all afterglows (from 
100 s to 100 ks), /ag = -Eend/£'start, IS auticorrelatcd with the 
isotropic-equivalent GRB output E^, for any model variant 
that can account for a significant fraction of the afterglow 
breaks listed in Table 3. 

The — /ag anticorrelation is illustrated in Figure 1 



* For /3x > 1/2, the fastest decay Ox = 2/3x + 3/2, exceeds the 
Ox = /9x + 2 limit for the large-angle emission decay, resulting 
when the emission from a spherical outflow stops suddenly; how- 
ever, for a jet seen after the jet-break time, the large-angle emis- 
sion does not exist because there is no radiating fiuid at angles 
larger than P"^ 
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Figure 2. Similar to Fig 1, but for energy-injection break as 
the origin of tlie X-ray liglit-curve break, with the afterglow out- 
flow being a conical jet and the jet-break having occurred at early 
times (before the first X-ray measurement). The ambient medium 
is homogeneous and X-ray assumed to be below the cooling fre- 
quency. Black points are for synchrotron as the dominant X-ray 
process, red points are for an inverse-Compton model. 



for two jet-break models and in Figure 2 for two energy- 
injection break models. For the synchrotron jet-break model 
(black points in Figure 1), the linear correlation coefficient 
in log-log space is r = — 0.42±0.08, yielding a probability for 
a stronger anticorrelation prob < 0.080 in the null hypoth- 
esis. For the inverse-Compton jet-break model (red points 
in Figure 1), r = -0.40 ± 0.07 and prob < 0.076. For the 
synchrotron energy-injection break model (black points in 
Figure 2), r = -0.54 ± 0.04, and prob < 0.52 x 10"^, while 
for the inverse-Compton energy-injection break model (red 
points in Figure 2), r = -0.43±0.06, and pro6 < 0.70 xlO"^ 
Similarly high chance probabilities are obtained for other 
models; thus the iS-, — /ag anticorrelation is only "tentative". 

The above linear correlation coefficients (Pearson's) r 
have been calculated assuming a 50 percent uncertainty (in 
natural scale) for both quantities plotted, the chance prob- 
abilities of those ant icorrelat ions have been estimated con- 
servatively for a linear correlation coefficient |r| — cr(r). The 
best-fits shown in Figures 1 and 2 were obtained by mini- 
mizing using the uncertainties of both quantities. 

That i) the fractional increase of the forward-shock en- 
ergy over a fixed time interval is determined by the expo- 
nent e of the power-law energy injection which accommo- 
dates the measured X-ray flux decay index Qx, and that ii) 
da/de < (a stronger energy injection yields a slower flux 
decay) suggest that the — /ag correlation may originate 
from a log — correlation. Indeed, as illustrated in Fig- 
ure 3, more energetic bursts are followed by faster decaying 
X-ray afterglows. The correlation between these two observ- 
ables is slightly stronger than the E^y — f^^ correlation for 
afterglows with jet-breaks at the X-ray light-curve epoch, 
but significantly weaker if the X-ray light-curve breaks are 
energy-injection breaks. 



The E., — /ag correlation implies that either the effi- 
ciency of the GRB mechanism or the energy of the ejecta 
producing the GRB emission are ant icorrelat ed with the ki- 
netic energy of the ejecta interacting with the forward-shock 
after the burst phase. 

If the dissipation mechanism leading to the produc- 
tion of 7-rays is internal shocks (Rees & Meszaros 1994) 
in a fluctuating relativistic wind, then the burst dissipation 
efficiency is determined by the ratio between the Lorentz 
factors of interacting shells. For an impulsive ejecta re- 
lease (short-lived engine), all the ejecta are produced quasi- 
instantaneously and the energy injection occurs as the in- 
ner parts of the outflow catch-up with the ever-decelerating 
forward-shock (Rees & Meszaros 1998). In this model, the 
E^ — /ag correlation implies an anticorrelation between the 
average shells' Lorentz factor ratio with the energy of the 
shells behind the leading edge of the outflow, smoother out- 
flows (i.e. with a smaller Lorentz factor contrast) having 
more energy localized toward the trailing edge. For an ex- 
tended ejecta release (long-lived engine), where the ejecta 
are released over a source-frame duration comparable to 
the observer-frame time during which energy injection takes 
place, the E^ — /ag correlation would imply an anticorrela- 
tion among bursts between the Lorentz factor contrast in 
the leading, GRB-producing part of the outflow and the ki- 
netic energy of the outflow expelled by the central engine at 
later times. This is a puzzling implication, which makes the 
long-lived engine scenario less likely. 

Irrespective of the dissipation mechanism, if the GRB 
output is correlated with the energy of the GRB-producing 
ejecta, then the E^ — /ag correlation implies a relationship 
between the GRB ejecta energy and the radial distribution 
of the afterglow outflow energy. For a short-lived central en- 
gine, where the GRB and afterglow outflows are the same, 
the Ery — /ag correlation would require that more energetic 
outflows have less energy behind the leading front. For a 
long-lived central engine, where the GRB outflow is the early 
afterglow blast- wave, the E.y — /ag correlation implies that 
engines that are initially more energetic release less energy 
after the burst. Perhaps this scenario appears more plausible 
because it opens the possibility of a constant total ejecta en- 
ergy. Unfortunately, we do not find that fa^gE-y is universal; 
instead, its distribution extends over at least two decades, 
for the /ag calculated in any light-curve break model. 

We note that the detection of X-ray fiares in some Swift 
afterglows indicates an extended ejecta release, but the en- 
ergetic output of those flares is only about 10 percent of the 
prompt GRB output (Falcone et al 2007). In contrast, the 
energy injection required during the afterglow phase corre- 
sponds to an increase of the outflow energy after the burst 
by a factor up to 10* (Figures 1 and 2). To reconcile the 
flare and afterglow energetics, the internal shocks that pro- 
duce fiares must be very inefficient X-ray radiators compared 
to the forward-shock (see also Maxham & Zhang 2009) or, 
else, the X-ray flares could be produced by the central en- 
gine, bearing little correlation with the energy of the ejecta 
released after the burst. 
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Figure 3. Correlation of power-law index Oxl of the early af- 
terglow X-ray flux decay and GRB output E-^ for 44 Swift af- 
terglows. Linear correlation coefficient is r = 0.34 ± 0.05, with a 
chance probability prob < 0.026 . 



4.2 Afterglow energetics 

The large increase of the forward-shock energy, by a factor 
f^g = 10 - 300 (Fig 1) or /ag = 10 - 10'' (Fig 2), required to 
account for the hght-curvc breaks of Swift X-ray afterglows 
raises the question if the total outflow energy is not too 
large and incompatible with the expected energy reservoir 
for a few solar-masses black-hole and its sub solar-mass torus 
formed in the collapse of a WR massive star. 

The jet energy can be estimated by assuming that, at 
the first X-ray measurement (tstart), the isotropic-equivalent 
forward-shock energy is the GRB output E^,. This assump- 
tion introduces an underestimation of the jet energy by a 
factor up to a few because energy could be injected into 
the forward-shock even before tstart (100 - 1000 s), maybe 
starting as early as the burst end (10 - 100 s), at a rate 
E oc t''^'^ — (this is the range of energy injection laws that 
accommodate the early X-ray flux decays of Table 3), and 
with the estimated jet energy satisfying Ej oc [E{tsta,rt)]^^'*' 
for a homogeneous medium and Ej oc [J?(tstart)]^^^ for a 
wind. For a short-lived engine, where all the outflow pro- 
duces the burst emission and only its leading edge yields 
the early afterglow, we have i5(tstart) ^ E^; in this case, 
the total jet energies calculated here may be substantially 
overestimated. For a long-lived engine, where the leading 
part of the outflow produces both the burst and early after- 
glow emissions, i5(tstart) = is equivalent to assuming a 
(reasonable) 50 percent GRB efficiency. 

The jet initial aperture 6j, needed to calculate the jet 
energy from isotropic equivalents, can be constrained as fol- 
lowing. For the model of a narrow jet whose edge becomes 
visible before the first X-ray measurement and the X-ray 
light-curve break is an energy injection, Oj < l/r(fstart). For 
the model where the X-ray light-curve break (at tb) is a jet- 
break undergoing a steady energy injection, 6j — l/r(tb). 
For the model of a wide jet whose edge becomes visible af- 
ter the last X-ray measurement (tend), > l/r(tend). Prior 



to the jet-break, the jet Lorentz factor is given by the dy- 
namics of a spherical blast-wave: 



r{t) = 800 



E53{t) 



no 



1/8 



t \ -3/8 



(27) 



for a homogeneous medium of particle density no per cm ^ 
and 



— r{t) = 175 



1/4 



t \ -1/4 



(28) 



for a wind-like medium of density n oc Ar normalized to 
that corresponding to progenitor mass-loss rate-to-terminal 
speed ratio of {dM/dt)/v = {Mq yr''^) / {W^ kms-^), 
respectively, with E-,[i{t) = 10"^"^ E{t) being the blast-wave 
isotropic-equivalent energy at observer time t. The latter is 
calculated assuming that i?(tstart) = E^ and for the energy 
injection law E act" required for each model to account for 
the measured X-ray flux decay. The observer time in the 
above equations corresponds to the arrival-time of photons 
emitted by the fluid moving at an angle 9 ^T"^ relative to 
the direction toward the observer, from where arises most of 
the high-energy emission received by the observer at time t. 

Putting together all the elements described above, we 
have calculated the jet final energy (at tend) for the X-ray 
afterglows whose light-curve breaks can be explained by the 
a jet-break model with steady energy injection or an energy- 
injection break model (the models listed in the mid and 
lower parts of Table 4). The jet energies for rio = 1 and A,, = 
1, averaged over all the afterglows that can be accounted for 
by each model (fewer than indicated in Table 4, because 
the jet energy calculation requires the GRB output/burst 
redshift), are listed in Table 5. For jet-breaks models, the jet 
opening is determined by the light-curve break epoch. For 
energy injection-break models, the epoch of the first (last) 
X-ray measurement provide an upper (lower) limit on the 
jet-break epoch, which yields an upper (lower) limit on the 
jet opening and, consequently, on the jet energy. Numbers in 
round brackets give the dispersion of individual jet energies 
(or their lower/upper limits). 

A potentially important result shown in Table 5 is 
that models where X-ray is inverse-Compton, the ambient 
medium is a wind, and X-ray is lower than the cooling fre- 
quency, can be ruled out on energetic grounds because, for 
them, the required jet total energy exceeds lo'''' " erg^. 
Those models lead to a high final jet energy because a sub- 
stantial increase of the jet energy is required to account for 
the measured X-ray flux decay rate, which stems from that 
the corresponding adiabatic jet model yields a steep X-ray 
flux decay (Table 2: = p+1/2 = 2/3x -|- 3/2, for a conical 



t Magneto-hydrodynamic simulations of accretion disks (e.g. 
McKinncy 2005, Hawlcy & Krolik 2006) show that the energy 
in relativistic matter (electromagnetic outflow) extracted from 
an accreting black-hole increases from 0.2 (0.03) percent of the 
accreted mass, for a non-rotating black-hole, to 15-40 (5-20) per- 
cent for a maximally rotating one. The maximal energy budget 
quoted here corresponds to a 20 percent MHD efficiency and a 
1 Mq accretion disk. 
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Table 5. Average (and dispersion) of the logarithm of jet energy 
(in ergs) at the last X-ray measurement, for the models in the 
two lower sections of Table 4. 



MODEL 


HOMOG MEDIUM 


Wind-like medium 




U < Uc 




1/ < Vc 




jet-break 


J-syn-cjet 


50.4(0.5) 


49.9(0.6) 


51.1(0.6) 


49.9(0.3) 


J-syn-sjet 


50.3(0.8) 


50.1(0.7) 


51.5(0.6) 


50.2(0.4) 


J-iC-cjct 


•50.7(0.6) 


49.8(0.7) 




50.0(0.3) 


J-iC-sjcl 


50. ()((). 7) 


l<).(i(().(i) 




19.7(0.3) 


energy-injection breaJc 


El-syn-cjet < 


50.0(0.6) 


49.7(0.6) 


51.7(0.7) 


49.9(0.5) 


El-syn-sjet < 


50.8(0.7) 


50.2(0.7) 


51.7(0.6) 


50.3(0.4) 


El-syn-sph > 


51.0(0.5) 


50.6(0.4) 


51.8(0.6) 


50.6(0.3) 


El-iC-cjct < 


50.2(0.7) 


49.7(0.7) 


53.6(0.8) 


50.0(0.5) 


El-iC-sjct < 


50.4(0.7) 


48.9(0.8) 


55.0(1.3) 


49.5(0.6) 


El-iC-sph > 


51.4(0.6) 


50.7(0.5) 


53.3(0.8) 


50.6(0.3) 



jet), from that energy injection has a weak effect in miti- 
gating the flux decay {da^/de = p — 2 = 2/3x — 1, for the 
spreading jet model), or a combination of those two factors 
(for the spherical outflow model). 

Assuming the same jet energy at the beginning of the 
afterglow, one can also calculate the ratios of jet energies 
among the three types of X-ray break models and outflow 
dynamics: i) a narrow-jet model with an energy-injection 
break (models E-cjet and E-sjet in Table 5), ii) a jet-break 
model with steady energy injection (models J-cjet and J- 
sjet), and Hi) a widc-jet model with energy- injection break 
(model E-sph). The interesting result is that, going from 
type i) to type ii) models above increases on average (i.e over 
afterglows) the total jet energy by a factor 2-3, with another 
increase by a factor 2-3 required when going from type ii) 
to Hi) models. These results hold for either type of ambient 
medium, either location of the cooling frequency relative 
to the X-ray, and for either X-ray radiative process, if the 
jet does not spread sideways, and only for inverse-Compton 
emission, if the jet is spreading laterally. Thus, the most 
economical model for a light-curve break is that of a narrow 
jet with an early jet-break and a broken power-law energy 
injection origin for the X-ray light-curve break. However, 
these conclusions do not hold for spreading jets if the X- 
ray radiation process is synchrotron. In this case, we flnd 
that the required jet energy is about the same for all model 
types above. Note that these relative energy requirements 
cannot be deduced from the average jet energies listed in 
Table 5. Here, we are referring to averages over afterglows 
of jet energy ratios between two models that explain the 
same afterglow break, while Table 5 gives jet energies for 
a given model averaged over all the afterglows that can be 
explained by that model. 



4.3 Loose ends 

The first caveat/ weakness of the model testing presented 
above was already discussed: the energy-injection break ori- 
gin for X-ray light-curve breaks is only weakly falsiflable 



by a comparison between the model a — ,5 closure relation 
and the observed flux power-law decay index a and spectral 
slope j3. The only way for this model to fail is an X-ray flux 
that decays steeper than expected for an adiabatic shock, as 
energy injection can only mitigate the afterglow flux decay. 

Furthermore, a comparison of the model light-curve 
breaks against observations is needed to validate possible 
X-ray light-curve break models, as the sharpness of X-ray 
light-curve breaks varies among afterglows. The majority of 
the well-sampled light-curve breaks used here are sharp, last- 
ing less than a factor two in time. Quite likely, the sharpest 
breaks cannot be accounted by a wind-like medium, what- 
ever is the break mechanism, owing to the slower jet deceler- 
ation produced by a decreasing external density (Kumar & 
Panaitescu 2000) . Without taking into account such limita- 
tions, testing only the ability of some models to account for 
the observed flux power-law decay indices will overestimate 
the fraction of afterglows that can be explained by a given 
model. 

We emphasize the importance of using afterglows whose 
light-curves have been sampled over a sufficiently long time 
before and after the break, so that the asymptotic flux de- 
cay indices can be measured accurately. Otherwise, shortly- 
monitored light-curves with breaks can lead to an underes- 
timation of the break magnitude 5a — a^2 — Qxi , which will 
favour a steady energy injection interpretation, because its 
effect is to reduce the break magnitude 6a. We note that 
eighty percent of the power-law segments used in this work 
last a decade or more in time. 



4-3.1 Optical afterglows 

Inclusion of the optical light-curves could provide a further 
test to the origin of X-ray light-curve breaks (e.g. Liang et al 
2007; Curran et al 2009; Liang et al 2008) in the two models 
discussed above (jet-break and energy- inject ion break). Both 
models yield achromatic light-curve breaks, appearing at the 
same time in all light-curves. For those afterglows where it 
can be established that the optical and X-ray emissions arise 
from the same source, either because the light-curve breaks 
are achromatic or because the optical and X-ray spectra are 
on the same power-law continuum, testing the light-curve 
break origin and determining the model details (location of 
cooling frequency, ambient medium stratification, radiation 
process, outflow dynamics) can be done by comparing the 
optical flux decay indices Qq with the model expectations, 
for the energy injection parameter e required by the X-ray 
flux decay indices ax. 

In GRBlog database of optical light-curves (Quiniby, 
McMahon & Jeremy 2004), we find five afterglows (050401, 
061121, 070420, 090424, 090618) listed in Table 3 with chro- 
matic X-ray breaks, hence their X-ray and optical emis- 
sions arise from different parts of the outflow, and eight 
with achromatic breaks, thus a comprehensive test of the 
X-ray break models cannot be performed at this time. We 
note that, of those 8 afterglows with achromatic breaks, two 
(060729, 081029) have Qx ^ Oo, for five (050922C, 051109, 
060614, 080710, 090424) 0.25 < Ox - Oo < 0.5, and for one 
(090510) ax - ao ~ 1. 

In general, ax = Oo (as for 060729) is expected if the 
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cooling frequency is not between the optical and X-ray and 
Qx — cvo = — (l/2)(ilogi'c/dlogi if Vc is in between. Accord- 
ing to Table 1, a decreasing Vc (which leads to ax > Oio), 
is obtained only for a homogeneous medium, thus all six 
afterglows with Ox — Oo > 0.25 require such an ambient 
medium. Energy injection accelerates the evolution of fc 
and increases the index difference ax — Oo above the ex- 
pectation for am adiabatic outflow: a^ — Uo = (e -I- l)/4 for a 
spherical/conical outflow and ao = e/?> for a spreading 
jet. Thus, the energy injection indices e in the range (0.3, 1) 
required by the X-ray decays of those six afterglows imply 
that 0.3 < ax — ao < 0, 5, consistent with what is mea^ 
sured for five of thorn, while the large fix — (lo measured for 
090510 indicates that tfie optical and X-ray are from differ- 
ent radiative processes (synchrotron and inverse-Compton, 
respectively) . 



4-3.2 Reverse shock 

In this work, we have considered only the forward-shock 
emission and ignored the emission from the reverse-shock, 
which energizes the GRB/post-GRB ejecta, the agent of en- 
ergy injection into the blast-wave. The two shocks dissipate 
equally the kinetic energy of the ejecta, which indicates that 
the reverse-shock could make a comparable contribution to 
the afterglow emission. 

For a short-lived ejecta source, it can be shown from 
the kinematics of the catch-up between the GRB ejecta 
and the decelerating forward-shock that the ratio of the 
Lorentz factor Fi of the ejecta arriving at the forward-shock 
at time t to that of the swept-up ambient medium T(t) is 
Fi/F(t) = [(4 - s)/(l -I- e)]^/^ ^ 1 (evidently, the ratio of 
those two Lorentz factors must be above unity, but cannot 
be much above unity, else the catch-up would have occurred 
earlier). In this case, the reverse shock is only mildly rela- 
tivistic and the ejecta electrons accelerated by the reverse- 
shock could be much less energetic than the electrons behind 
the forward-shock, making the emission of the latter domi- 
nant at X-ray photon energies. 

For a long-lived central source, the contrast between 
the Lorentz factor of the unshocked outflow and that of 

the forward-shock can be arbitrarily large, the reverse-shock 
could be highly relativistic, provided that the ejecta are 
not too dense (e.g. Sari & Piran 1995). The decay of the 
reverse-shock flux can be calculated in a way similar to that 
shown here for the forward-shock; however, the parametriza- 
tion of the reverse-shock emission is not that simple because 
i) its dynamics depends on the Lorentz factor of the ejecta 
and their density, and ii) its peak flux depends on the total 
ejecta mass having entered the reverse-shock. Thus, the cal- 
culation of the reverse-shock emission requires at least two 
parameters that quantify the energy and mass mass injec- 
tion into the blast-wave, leading to a model even less con- 
strainable (observationally) than the forward-shock (which 
requires only one parameter, e, for the energy injection law). 



5 CONCLUSIONS 

By comparing the flux decay index (a) with the spectral 
slope {P) measured for a set of 98 well-sampled Swift X- 
ray afterglows with the expectations for the forward-shock 
model, we find that about a third of those afterglows display 
a break that could be the traditional 'jet-break' claimed to 
have been observed in a dozen pre-Swift optical afterglows 
(§2.3.1). For this test of adiabatic jet-breaks, we have con- 
sidered various details of the forward-shock model: homoge- 
neous or wind-like ambient media, spectral cooling frequency 
higher or lower than the X-rays, synchrotron or inverse- 
Compton as the dominant X-ray process, a coUimated con- 
ical or spreading outflow. 

If a steady energy injection into the forward-shock oc- 
curs, then the jet-break model may account for 88 percent of 
the 117 light-curve breaks or, else put, for at least one break 
in 87 percent of the 98 afterglows selected here (§3.3.1). This 
vast improvement over the adiabatic jet model suggests that 
energy injection could be a prevalent process in GRB after- 
glows. 

The remaining 12 percent of breaks that cannot be ex- 
plained by a jet-break undergoing a steady energy injec- 
tion, and the 19 afterglows with two breaks (only one be- 
ing, at most, a jet-break), indicate that some X-ray breaks 
may originate from a change in the rate at which energy 
is added to the forward-shock (an energy-injection break), 
and it is possible that some (many ?) of the identified jet- 
breaks are also energy-injection breaks (§3.3.2). It should be 
recognized that the complete success of the energy-injection 
break model over the traditional jet-break model could be 
due mostly to that the latter model is over-constrained (ob- 
servations provide two constraints: flux pre- and post-break 
power-law decay indices for one model parameter: the en- 
ergy injection power- law exponent), while the former model 
is not. 

We find a tentative correlation between the GRB out- 
put and the afterglow X-ray flux decay rate, more energetic 
bursts being followed by faster decaying afterglows (§4.1). 
If energy injection is a common process in GRB afterglows, 
then the slower decaying afterglows should be identifled with 
a stronger total energy injection. This means that the GRB 
energy - afterglow (X-ray flux) decay rate correlation should 
lead to an anticorrelation of the GRB output with the in- 
crease of the outflow kinetic energy during the afterglow 
phase, i.e. more energetic bursts should be followed by af- 
terglows with less energy being injected. This anticorrelation 
is, indeed, observed for both the jet-break and the energy 
injection break models for X-ray light-curve breaks, and for 
various model details (dynamics, X-ray afterglow emission 
process, radial stratiflcation of the ambient medium). If the 
GRB progenitor is a long-lived source of relativistic ejecta, 
the above anticorrelation implies an anticorrelation of the 
kinetic energy of the leading outflow (producing the burst) 
with that of the following ejecta (which inject energy in the 
blast- wave during the afterglow phase), but the total ejecta 
energy (burst plus afterglow) is far from being universal. 

The energy budget required by most of the possible vari- 
ants of the forward-shock model ranges from 10*® to 10®^ 

ergs, and only variants involving inverse-Compton emission 
in the X-rays and a wind-like medium can be excluded on 
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energetic grounds because, in those models, the required jet- 
energy exceeds 10^^ erg (§4.2). Considering both models for 
light-curve breaks (jet-break with steady energy-injection 
and energy- injection break), we find that the most ener- 
getically economical model is that of a narrow jet whose 
jet-break occurs early on (before first X-ray measurement) 
with the light-curve break being an energy-injection break, 
followed by a model where the X-ray light-curve break is 
a jet-break, with the most "wasteful" model being that of 
a wide jet whose break occurs late (after last X-ray mea- 
surement), the X-ray break being an energy-injection break. 
However, the average (over afterglows) energy ratio among 
pairs of these three break models is only 2-3. 
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